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Abstract 

Hylleraas-Configuration Interaction (Hy-CI) calculations on the ground 1 S state of helium atom 
are presented using s-, p-, and d-Slater orbitals of both real and complex type. Techniques of 
construction of adapted configurations, optimization of the orbital exponents, structure of the 
wave function expansion are explored. A new method to evaluate the two-electron kinetic energy 
integrals occurring in the Hy-CI method is presented. The calculations show that nanohartree 
accuracy, about 0.0002cm -1 is achieved. 
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1. Introduction 



The Configuration Interaction method (CI) is of great importance in quantum mechanical 
calculations of the electron structure of atoms and molecules. It is well known that the 
shortcomings of the CI method are due to the form of the wave function which does not 
fulfill the electronic cusp condition 
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The CI wave function does not contain explicitly odd powers of the interelectronic coordinate 
Tij (implicitly the CI wave function does contain terms rf •, rf -, • ■ ■ rf n formed by combination 



nn n 

of angular orbitals p, d, f |2|-|4|). But odd powers are the ones energetically important [5[. 
In the CI wave function the piling up of higher and higher angular terms attempting to 



represent the term like in a Tailor expansion [5| shows how important the details of 
the cusp are. The energy improvement when increasing the quantum number / follows an 
asymptotic formula proportional to f i + 1/2)- Q in the case of he.ium atom, and in genera, 
for a larger number of electrons p, |8J. The nuclear cusp condition is always fulfill, in the CI 
as in the Hartree-Fock (HF) wave functions: 

~) = ~Z, (2) 

Z is the atomic charge, or the orbital exponent. The cusps (positive for repulsion and 
negative for attraction) account for two-body correlation, but not for three-body correlation. 
These conditions are a result of the singularities of the Hamiltonian at fj = and = 0: 

H _ lyd 2 y}_d__yl + J L 2 d I y rf + r 2 12 - Q 2 
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As the exact wave function is obtained from the equation: = E. This equation 

leads to the exact energy only if the cusps conditions of Eqs. (1,2) are fulfilled. The CI 
wave function is then not a formal solution of the Schrodinger equation and therefore the 
CI wave function converges very slowly. Schwartz jt| has demonstrated that the inclusion 
of other tern^thaa r 12 into the wave function, like negative power, of r, fl, fractionary 
powers of [11( and logarithmic terms ln(r, + rj) accelerates the convergence of the wave 
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function to the exact solution. Logarithmic terms are important to describe the three- 



particle coalescence region 
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is] ]. Nakatsuji [l^] has proposed the Iterative-Complement- 
Interaction (ICI) method which generates those terms which are necessary in the structure of 
the wave function for very highly accurate calculations. Including a ln(r x + r 2 + /3r 12 ) in the 
wave function a 40 decimal digits accurate energy was calculated, showing that logarithmic 
terms are necessary. Such accurate calculations are restricted to the case of two-electron 
systems. 

Other of the shortcomings of the CI wave function are the functions used as orbitals. 
If the functions form a complete set it is obvious that the exact solution can be expressed 
in terms these functions. Some types of orbitals form a complete set (or overcomplete in 
case of double basis) like Slater orbitals and Sturmians. If the set is not complete the wave 
function needs some states from the continuum. The importance the continuum functions 
in representing the ground state of helium atom has been pointed out by many authors 
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16j. If one carries out calculations to the ultimative limit it is immaterial which set is 



used. Therefore an infinite expansion of CI configurations would be the exact wave function. 
The reason why the CI wave function gives good results is the consecutive inclusion of angular 
functions which represent r^. 

It is then clear that the explicit correlated wave functions including the interelectronic 

n 

coordinate proposed by Hylleraas [17J are an alternative to the CI wave function. Sims 
and Hagstrom 18|, [19| introduced the Hylleraas-Configuration Interaction (Hy-CI) wave 
function which combines the use of higher angular momentum orbitals as CI does with the 
inclusion of the interelectronic distance into the wave function as the Hylleraas-type wave 
functions. The first terms of a Hylleraas-CI wave function are CI terms. 

Helium atom has been subject of numerous investigations. Its electronic structure has 
been determined to the greatest accuracy known in quantum chemistry. Our interest is to 
test in the case of helium atom our integral subroutines and computational techniques as 
the spin, antisymmetrization, construction of symmetry adapted configurations, the usage 
of the Hamiltonian in Hylleraas coordinates, and to win experience to investigate larger 
systems. 

This work will be thoroughly referred and compared first, to the early helium calculation 
of Weiss |5| for which there is enough computational data, and second, to the highly accurate 
more modern Configuration Interaction and Hylleraas-Configuration Interaction (Hy-CI) 
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calculations of Sims and Hagstrom j^]. The CI calculation of Weiss is taken as starting 
point to check our program code. The CI and Hy-CI calculations of Sims and Hagstrom 



have been taken for comparison as t 



calculations of Sims and Hagstrom 



re s,p,d-limits. Finally we have reproduced Hy-CI 



20). 



Novel in this work is the use of a different technique to evaluate the kinetic energy Hy- 
CI integrals which we have used to reproduce the energy values obtained from Sims and 
Hagstrom. 

2. The Hylleraas-CI wave function 



The Hy-CI wave function 19| for the 1 S ground state of helium atom is written: 



N 

V = Y, C k®k, $k = 6(L 2 )A^ kX , X = M - Pa), (4) 

k=l 

where are symmetry adapted configurations, N is the number of configurations and the 
constants Cj~ are determined variationally. The operator 0(L 2 ) projects over the proper 
spatial space, so that every configuration is eigenfunction of L 2 . A is the 2-particle anti- 
symmetrization operator, and x is t ne spin eigenfunction. The Hartree products ipk, are 
products of orbitals of the required symmetry multiplied by the interelectronic coordinate: 

ipk = 0jfc(l)0ifc(2)ri 2 , v = 0, 1 (5) 
For v = the Hylleraas-CI wave function reduces to the CI wave function. The power 



22J, we know that 



v = 1 does not represent any restriction. As discussed in Refs. 
even powers of the interelectronic coordinate are equivalent to products of p-, d-, ■ ■ ■ type 
orbitals, for instance: 
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= p(l)p(2). (6) 



Furthermore, higher odd powers can be expresses as r 12 r^- It can be demonstrated 21] 
that 7*12 is an infinite expansion of angular orbitals: 

s(l)5(2)r 12 = a(l)8(2) + p(l)p(2) + d(l)d(2) + /(l)/(2) + . . . (7) 

In the case of two-electrons systems Sims and Hagstrom [3 1 ] have shown that the Hy and 
Hy-CI wave functions are equivalent. 
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3. The Hamiltonian 



For a Hy-CI wave function in helium atom the effective Hamiltonian in Hylleraas coor- 



dinates may be written [35] 



23|: 
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Using the Hamiltonian in this form is easy to derive the involved integrals. The angular 
momentum operator can be extracted and its eigenvalue equation used: 



d 2 1 A cot ft d 



2j^r 2 d6 2 2 f^r 2 sin 2 d i d(p 2 2 rf dO- 



(9) 



L 2 4>i = + 

with li the angular quantum number of the orbital </>;. 

From the variational principle one obtains the matrix eigenvalue problem: 



(10) 



where the matrix elements are: 



fH — EA)C = 0, 



H kl = / QkH^idr, 



(11) 



A 



hi 



(12) 
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4. The use of real Slater orbitals 

We have constructed a set of real s-, p-, <i-Slater orbitals which are orthogonal and 
unnormalized. The exponents are considered as adjustable parameters. They are defined 
as: 



ns = r n 1 e ar 



np z = r n_1 e _/3r cos(#) 
np x = r n ~ 1 e~ l3r sin(0) cos(</?) 
np y = r n_1 e _/3r sin(6 l ) sin(<^) 
nd 2 2 = r n - 1 e - 7r (-cos 2 (9) ' 



2 2, 
nd xz = r n ~ l e~ ir sin(6 l ) cos(#) cos((p) 

ndy Z = r n_1 e~ 7r sin(0) cos(6 l ) sm(cp) 

nd xy = r n ~ 1 e~ 1T sin 2 (#) cos(y>) sin(yj) 

nd x 2_ y 2 = r n - 1 e- 7r sin 2 (^)(cos 2 (^) - sin 2 (</?)) (13) 
The orbitals are eigenfunctions of L 2 , but they are not eigenfunctions of L z . 

4.1 CI Integrals over real Slater orbitals 

The evaluation of the resulting matrix elements Eq. (11) for a CI wave function leads to 
)roducts of one- and two-electron integrals. These integrals can be evaluated as shown in 



24j , we write here compact expression of all non- vanishing two-electron integrals 36[ over 



p-orbitals: 
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(8(1)8(2)^-8(1)8(2)) = (47T) 2 / — 

( S (l) Pi (2)— 8(1)^(2)) = ^/— ). / : .r./,.: 
ri2 3 \ ^12 , 



(Pi(l)a(2)— a(l)ft(2)) = ), • = x,y, 2 

ri2 9 \ ^ 2 7 



(^(1^(2)^(1)^(2)) = M 2 /^ - * ^ •?> ^ = 

( ft (l) Pj (2)i-p,(l) Pi (2)) = (4vr) 2 ^ J V i ^ j, 



i,j = x,y,z (14) 



where S12 is the smallest of r\ and r 2 , and <?i2 the largest of r± and r 2 . The radial part of 
these integrals is expanded in terms of auxiliary integrals 24{. The non- vanishing integrals 
including <i-orbitals are: 
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(a(l)du(2)— s(l)d iz (2)) = ^-( S ~f), i = z,y (15) 
r 12 to \gi2/ 

[s{l)d zz {2)—s{l)d zz {2)) = MY J-\ (16) 



5"12 25 \5 , 12/ 

(a(l)du(2)— du(l) a (2)> = V i = x,y (17) 

r i2 to \g 12 / 

[8{l)d„(2)±d„(l)s(2)) = ^/#Y ( 18 ) 



r 12 25 \ gf 2 

(B(D^(2)^(1)^(2)) = (4ir)»(l-L + ), i = *,y (19) 

(^(1)^(2)^(1)4,(2)) = (4.) 2 (li- - ^|). (20) 

( ft (l)^(2)i-^(l) Pi (2)) = (4vr) 2 ^g + 3^g^), * = *>V (21) 

^(1)^(2)^(1)^(2)) = (^) 2 (^| - (22) 



^(1)^(2)— p,(l)d„(2)> = ±-L( -f ), ^ = x,y (23) 



(47r)_ , 

(p ,( 1)(i „(2)i- d „(l)p I (2)> = (4.) 2 (-^| + I ^|), i = *,» (24) 



1 _ „w ,„» (4 



2 / „2 
^12 



(p x (l)d xz (2)— Py(l)d yz (2)) = ^-(J? ), ( 25 ) 
1 , ™ (4vr) 2 / s ?2 



(p x (l)d xz (2)—d yz (l)p y (2)) = '—t-< -f ), (26) 

7"i2 MO \<7 12 , 



(^(1) S (2)— pi(l)p a (2)) = ^/^r ). » = X >V (27) 
n 2 45 \gf 2 ' 



(dUl) Pi (2)—s(l)p z (2)) = ^(4Y i = *>V (28) 
n 2 75 \^ 2 / 



(Pi(l)d iz (2)— p z (l)d«(2)) = ( ), z = x,y (29) 



(r(1) A (2)— ^(1)4,(2)) = ), i = x,y (30) 



^12 1 45 \# 2 2 /' 

These expressions have been obtained with the help of the algebraic computer program 
Maple |25j. 



S 



4.2 Construction of symmetry adapted configurations 

The ground state configuration of helium atom in the spectroscopic notation is of the 
type ss and has S symmetry. We can construct configurations using p-, d- orbitals of S 
symmetry. These configurations are: 



ip s = ns(l)ns(2) = ss = s 2 , (31) 
ip p = np{l)np{2) = pp = p 2 , (32) 
ip p = nd{l)nd(2) = dd = d 2 , (33) 



with 



i ) v = VV = Px(l)Px(2) + p y (l)p y (2) + p z (l)p z (2) = p x p x + p y py + p z p z . (34) 

The basis functions p y p y and p z p z may not contribute energetically as much as p x p x but 
they are necessary to have the proper symmetry, so that the wave function is eigenfunction 
of L 2 . Therefore a matrix element between two of these configurations has to be calculated 
as: 



pp 



H PP) = (p x p x + PyPy + P^ z \H\ P X P X + PyPy + ^ ^ Z ) . (35) 

One has to note that using symmetry adapted functions or configurations (SAF) and non- 
SAF the dimensions of the //-matrices are different, whereas a proof of the correctness of 
the calculation is that the diagonalization of the //-matrices in both cases should lead to 



the exactly same energy results 
or powers may be used, pp' . 
In the case of d-orbitals: 



20] . We have used the notation pp, but different exponents 



V'd — dd — d z 2d z 2 + d xz d xz + dy Z dy Z + d X yd X y + d x i_y2d x i_y2. (36) 
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5. The use of complex Slater orbitals 



The complex Slater orbitals with quantum numbers n, m and / are defined by an unnor- 
malized radial part and an angular orthonormal part which is a spherical harmonic: 



/(r) 



y.n-lg-arym*/ 



,n'-l„-a'rv"i' 



1 V 



(37) 



The spherical harmonics in Condon and Shortley phases 26|, p. 52] are given by: 

'21 + 1 (I -m)!l 1/2 



y; 



iT(cos0)e im * 



(38) 



4tt (/ + m)! 

where P z m (cos#) are the associated Legendre functions. The spherical harmonics and asso 
ciated Legendre functions used along this work are written explicitly in 0, p. 14]. 
The charge distributions are: 



i+v 



^L,M(r) = <f(r)0'(r) = ^ (2L + l) l ' 2 C L {l', m>; /, m^V^ 

L=\l-l'\ 



(39) 



where N = n + n' — 1, M = m' — m and the exponents u; = a + a'. N,L,M are the 
quantum numbers of the charge distributions. C L (l',m'; l,m) are the Condon and Shortley 



coefficients 



26|, p. 52]. 



5.1 Integrals over complex Slater orbitals 

The two-electron integrals appearing in the Hy-CI method are then defined: 

I(N 1 ,N 2 ;u 1 ,u 2 ;u)^ l 'r i l'^ 2 ' m2 = J fiv 1 ,L 1 ,M 1 (ri)fiv 2 ,L 2 ,M 2 (r 2 )ri2'irirfr2, (40) 
with the cases v = —1, 0, 1, 2. After angular integration the two-electron integrals are: 



I{N 1 ,N 2 ;u 1 ,u 2 ;u) lul , i / 2A 



mi ,m lt m2,Tn2 
i 

2 



-l) Jvfl 5(M 1 + M 2 ,0) ^ ^ 5(L 1; L 2 ) 

ii = Ki-«il L 2 =|J 2 -i^| 

2 

x J](2Z, + l) 1 ^ m<; ^ m,)!^, iV 2 ; Wl , w 2 ; 1/; L 2 ) (41) 

i=i 

where I(iVi, iV 2 ; wi, Co> 2 ; z/; L 2 ) are the basic radial two-electron integrals. The occurring basic 
two-electron integrals in the calculations of helium atom are: 
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I(Nx, N 2 - u u u 2 ; 0; L) = 6(L, 0)A(N 1 + 1, u x )A{N 2 + 1, u 2 ) 



(42) 



/(JVi,JV 2 ;wi,a; 2 ;l;L) 



1 



(2L + 1) 



x 



1 



( -> L _1) ^ V ( Nl +L + l,N 2 -L + 2;u 1 , u 2 ) + V(N 2 + L + 1, JV X - L + 2; w 2 , wi)] 
1 



+ 



(2L + 3) 
/(JVi.JVajWi.Wa;-!;^) 



[V(iVi + L + 3, N 2 - L; u h u 2 ) + V(N 2 + L + 3, Nt - L; u 2 , wi)] 
1 



(43) 



(2L + 1) 



[Vim + L + 1, N 2 — L; ui, oj 2 ) + V{N 2 + L + 1, N x — L; u 2 , 

(44) 



I(N U N 2 ; ut, io 2 ; 2; L) = 5(L, 0) [A(N t + 3, wi) A(iV 2 + l, u 2 ) + \A(N 2 + 3, w 2 )A(A^ 1 + 1, Wl )] 

- l)A(iV 1 + 2,o; 1 )A(iV 2 + 2, w 2 ) (45) 



A(n, a) and V^(/c, Z; a, (3) are auxiliary two-electron integrals denned in Ref. [29| . 



5.2 Hylleraas-CI two-electron kinetic energy integrals 



In this section we treat the two-electron kinetic energy inte gra ls appearing in the Hy-CI 



method when using the Hamiltonian in Hylleraas coordinates 



23] . For any atomic number 



N > 3 the kinetic energy integrals are of two- and three-e 



ectron type. The three-electron 



kinetic energy integrals have been evaluated in Refs. [281. |29I| . For N = 2 the integrals 
are of course, of two-electron type. Sims and Hagstrom [28| evaluated the kinetic energy 
integrals using the transformation of Kolos and Roothaan 30] which partially avoids the 
differentiation with respect to terms appearing on the right hand side of the matrix 
elements. In this work we will perform the derivatives directly over and using Hamiltonian 
written in polar and interelectronic coordinates, we shall prove its correctness and will solve 
the resulting integrals in terms of basic two-electron integrals defined in Ref. 29]. 

The Hamiltonian in Hylleraas coordinates can be separated in a sum of one electron 
potential and kinetic energy operators: 



H{1) = f (1) + V(l) = T R (1) + f e (l) + V(l). 



(46) 



i=i 
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The the kinetic energy operator can be again separated into radial and angular parts. Let 
us evaluate the kinetic energy of electron 1: 

f (1) = f Rl (l) + f R2 (l) + f Rs (l) + f L (l) + f dl (l) + T fc (l) + 4(1). (47) 
The radial operators are: 



and the angular ones: 



Id 2 Id 

T «^—2^--^' <48) 



TrM = (49) 
r 12 or 12 

rf m 1 r\ + r\ 2 - r\ d 2 

TrAI) = "2 nr 12 dndr^ (50) 



2 r{ 89 f sm^i <90f <9#i 

~ r 2 cosfl 2 9 2 

t\Tyi sint/i oOidr 12 

2 rfr i2 o9idr 12 

4(1) = --^^smCfc-fc)^^ , (54 ) 

the potential energy operator is: 

1/(1) = - A (55) 

as Z = 2 for helium. In two-electron kinetic energy integrals we have basis functions CI type 
with v — and basis functions Hy-CI type with v — 1. In the following, we will use z/ for 
power of r\2 in the left hand side basis functions and v' on the right ones. The interelectronic 



distance in the integral has then the power r 



v+v' 
12 • 



The nuclear attraction potential energy is: 



I PE (1) = (0(^)0^2)^2^10(^)0^2)^;) 

= -2/(iV 1 - 1, N 2 ; wx, wa; 1/ + ^ST™*' ^ = °' 1 ( 56 ) 
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with 2 the atomic nuclear charge. 

For the kinetic energy we have to evaluate: 

I KE (1) = (<K ri )0(r 2 K 2 |T R + f,|0( ri )0(r 2 )rr;) (57) 

The radial part: 

Ike,r = (H^)<f>(^y i2 \f Rl (l) + f R2 (l) + 7^(1)1^)^2)^) (58) 

performing the derivatives given by the operators Eq. (48-50) and using the definition of 
the two-electron integrals Eq. (40), it can be straightforward obtained: 



[n 



l2 



<-KE,Ri 



n l) t/i\t o a r 1 /\mi,m' lt m2,m' 2 

I(Nx - 2, N 2 ; wi, w 2 ; 1/ + v ) h ^J 2 



2 

+ (nWi + a'MNi - 1, N 2 ; u> u u> 2 ; v + v'^Xm 1 ^ 



/\mi,m 1 ,m2,m2 
1 
2 



a l T/ AT AT 1 i\m 1 ,m' 1 ,m 2 ,m' 2 

- ~Y I( y N ^ N ^ w i> ^ v + v ) hA / 2A \ 



i2 



(59) 



IKE,R 2 — 



v 
'2 



I{N u N 2] u u u 2 -v + v -2) hA ,il,v 2 

T/AT AT o 1 /\m 1 ,m' 1 ,m2,m' 2 

+ /(iVi, N 2 - 2; wi, w 2 ; v + 1/ ^ 

r/ A7" o AT 1 O 1 / o\ m l' m 'l' m 2> m 2 

- /(iVi - 2, iV 2 + 2; wi, u 2 - v + z/ - 2),^,,^,, 



+ 1/ 



2 L 



,mi,m 1 ,m2,m2 

1*2 



T / AT 1 1 AT 1 / o\ m l> m 1 1" 

/(iVi + 1, iV 2 ; u u u 2 \v + v - 2),^,,^,, 

r/Ar 1 at 1 /\mi,m!,m2,mi 

+ /(iV 1 -l,7V 2 ;a; 1 ,a; 2 ;i/ + i/) Ji)JliJ2 1 )J , 

—I{N\ -1,N 2 + 2; Wl , W2 ; 1/ + v' - 2)™$'™'^] , (60) 



Wi = -^/(JVi, iV 2 + 2; Wl , W2 ; 1/ + v' - 2)^'™'™*. (61) 

The factor v' ensures that CI configurations will not contribute and but the Hy-CI will 
contribute: v' — 1. The integral containing the the square of the angular momentum operator 
can be easily evaluated by means of its eigenvalue equation: 
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l KE,e u L 



:ri)0(r 2 )rr 2 |^|0( ri )0(r 2 )rr; 
2 r{ 



mi ,mj ,m2,rri2 



-I'Si + l)/(JVi - 2, iV 2 ; Wl , w 2 ; z/ + z/, L 2 )™J fa? 

/ i> i) ti 2 



(62) 



2 J- "> J- ' ' \ A ' «" «" ' '"'»l,»i,t2)»2 

where ^ is the quantum number of 0(ri). 

The evaluation of the angular kinetic energy is more involving. The following integrals 
should be evaluated: 



l KE,01 



^ri)0(r 2 )r 



r 2 cos# 2 <9 2 
riri 2 sin Q\ dd\dr\ 2 



(ri)0(r 2 )< 2 



(63) 



l KE,02 



( <v»2 ~,2 ~,2 

-J ( 0( ri )0(r 2 K 2 |^^ ^cotd 



5 2 



r?ri 2 



l KE,(f>3 



, \ if \ v I r 2 sin #2 

v ri)0(r 2 )r 12 | sm - </> 2 j 



(96»i <9r i 2 
<9 2 



^ri)0(r 2 )r 



V 

12 / 5 



(ri)0(r 2 )r 



(64) 



(65) 



rxr 12 sin 6*i d<pidri 2 
First the derivatives of the spherical harmonics with respect to the polar angle [31, 5.7.] is 

made: 



d6 1 



'l, vi, 



(66) 



and the cosine recursion relation is used: 



cos 9 2 Y,7' 



7 2, V2 



{l' 2 + m' 2 ){l' 2 -m' 2 ) 



1/2 



+ 



(2/ 2 + l)(2/ 2 -l) 
(/ 2 + m 2 + l)(/ 2 -m 2 + l^ 1/2 



(67) 



(2/ 2 + l)(2/ 2 + 3) 

then the products of spherical harmonics are linearized obtaining an expression in terms 



of Y, 



Mi, 



Li - l t7i,0i),y L2 - l f7 2 ,^ 2 



M 2 , 



I. After a rotation, see Ref. 



M 2 , 



L 2 WUVD 



29|, Eq. (22)] Y^ 2 (9 2 ,(/) 2 ) leads to 
Now the spherical harmonics with the same arguments are linearized again in 



order to avoid sigularities when using the recursion relations containing a sine function: 



n,(pi, 



E 

L=\L X -L 2 \ 



2L + 1 
Air 



1/2 



C (L 2 , M 2 ; Li, Mi)Y£ 



M, 



n,<pij 



(68) 
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C L (l' ,m';l,m) are the Condon and Shortley coefficients 26l. Egs. (6-11)]. The recursion 
relation containing the sine functions requests that M^O [29_|, 5.7.]: 



yM 



sin^i 



2M 



(2L + 1) 
(2L + 3) 



1/2 



[[(L - M + 2){L - M + l)] 1 /^^ 1 ^, 0x; 



+ ((L + M+l)(L + M + 2)) 1/2 e-^ i r i J| 1 (^ 1 ,0 1 )] , M+1,M-1^0. (69) 

this is fulfill in our case, because M = M x + M 2 + 1 or M = M 1 + M 2 - 1 and M 1 = -M 2 
due to the form of the symmetry adapted configurations in helium atom. Substituting into 
Eq. (63) we have: 



i[+h 

Li = \l[-h\ 



2L X + 1 



x 



{l' 2 + m' 2 ){l' 2 -m' 2 ) 



47T 
,1/2 



-| 1/2 ( l' 2 -l+h Li+L 2 

£ ( 2 )(2L 2 + 1) 1 / 2 £ ^ 

L 2 = K^-l-« 2 | L=|Li-L 2 | 



2L + 1 



2L + 3 

C La (Z 2 -IX; fe, m 2 )/(Zi, mi, Li, L 2 , L, M X) M 2 , M) 



1/2 



(2Z 2 -l)(2Z 2 + 3) 

x J(JVi -l,JV 2 + l;wi,w 2 ;v + ^-2;L2 

« 2 +l+i 2 ii+ij 
+ £ (2) (2^2 + l) 1/2 £ (2) 

L£=|^+l-i 2 | L'=|Li-L^ 

x C^+l, m 2 ; Z 2) m 2 )/(Zi, m' l5 L 1; L' 2 , L', M 1; M 2 , M)J(JVi -1, iV 2 + l; Wi, w 2 ; u+u'-2; L' 2 



~2L' + 1~ 


1/2 


r(Z 2 -m 2 + l)(Z 2 + m 2 + l)" 


1/2 


_2L' + 3_ 




I (2Z 2 + 3) 





(70) 



with the factor / is a function of the quantum numbers: 
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f{l' 1 ,m' 1 ,L l ,L 2 ,L,M u M 2 ,M) = 

(/■7T /"27T 
[(L-M l -M 2 + l)(L-M l -M 2 )] l/2 J J Y^\ +M2 (9 1 ^ 1 )sm9 1 d9 1 d ( p 1 

^ M2+2 (^i, ^Oe- 2 ^ 1 sinM0i#i) 
-(-l)^-i— ^— T i [(/; - mi + l)(Zi + m' x )] 1/2 C^'i, mi-1; J lf mOC^L., M 2 ; L 1; -M 1+ l) 

(PIT p27T 
[{L - M 1 - M 2 +3){L-M 1 - M 2 )] 1/2 J J Y^\ +M ^ 2 (9 1 , <\> x )e 2i ^ sin 9 x d9 x d<\> x 

Y^+ M ^e u 00 sin 9 1 d9 1 d(f )l \ . (71) 
The first kind of integrals over spherical harmonics are: 



7T p2ir 

M 1 +M 2 




0i) sinMM^i = 5(L + 1, 0)<J(M! + M 2 , 0). (72) 

'o 

The following integrals are evaluated integrating over the associate Legendre functions : 



7T p2tT 




K i A J 1 +M2+2 (0 1 ,(/. 1 )e- 2 ^ 1 sin OxdOxdfa = 2ir / P^+^+^cos^) shiM^i 

Jo 

n2ir pit 
Y^ 1 +M2 - 2 {9 l) (j) l )e 2i ^sixi9 l d9 l d(l) l = 2n P^ l +M2 ~ 2 {cos9 l )sixi9 l d9 l (73) 
Jo 

using the algorithm developed by Wong [33| . The overlap integral over Associate Legendre 
functions is: 



/•7T Pi max P2 max 

/ P^(cos9)P^(cos9)sm(9)d9=J2 E 



Pi P2 

a L- L M 1 a L 2 M 2 



Pl=0 P2=0 



r (|(Li ■ /,, - .U : - M, 2gi - 2p 2 + 1)) r (|(M! + M 2 + 2gi + 2p 2 + 2)) 

r(i(L 1 + L 2 + 3)) ' (74j 

T are Gamma functions. This formula is valid for < Mi < L\ and < M 2 < L 2 . If M 
is negative the formula is used: 

P^(cos^) = (-l) M §^|^ M (cos0i). (75) 
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For even L\ + L 2 — M\ — M 2 the integral vanish. The coefficients are: 



-1)p(L + M)\ 



(76) 



lm 2 m + 2 p{M + p)\p\{L - M - 2p)V 
To evaluate the second contribution to the angular kinetic energy we proceed as showed 
above and need the recursion relation [31, 5.7.]: 



-cot0Y^(0i,0i; 



2M 



\{L + M + 1)(L - M)] 1/2 e-^r L M+1 (0i, 0x) 



f— [(L - M + 1){L + M)] 1/2 e i ' t ' 1 Y^~ 1 (9i, fa). 



2M 



obtaining: 



(77) 



2M 



+ 



ft. 2 = ^ E (2) (2^i+l) 1/2 £ ^(2L 1+ l)V^(^,m-/ 2 ,m 2 ) £ ^ 

il=|«i-ii| i2 = Ka-'2| L=|Li-L 2 | L 

( _l)M 1+ il [(z / +m / + _ m /)]i/ac7^(/iX + 

n27T 
l^\ +M2+2 (0i, 0i)e- 2 ^ sin0id0i#i 

n27T 
1^+^(01, 0i) sin 

" (-l) Ml ~^[(^ -™i + 1)W + m' 1 )] 1 / 2 C Ll (/i,m' 1 - 1;^) 

/*7T /•27T 



,1/2 



2M 



+ 



^[(L + M x + M 2 )(L - Mj - M 2 + l)] 1 / 2 jT jf * ^+^(0!, 00 sin0 1 d0 1 # 1 

— [(L-M 1 -M 2 + 2)(L + M 1 + M 2 -1)} 1 / 2 J^ J Y^ M ^ 2 {e u fa)e 2i ^ sin 

x [J(JVi - 2, iV 2 ; W!, w 2 ; 1/ + 1/ - 2; L 2 ) - /(JVi, iV 2 ; w x , w 2 ; 1/ + 1/ - 2; L 2 ) 

- J(iVi - 2, iV 2 + 2; Wlj co 2 ] u + iS-2; L 2 )] , (78) 



The 



derivative of the spherical harmonic with m[ = with respect to <p\ vanishes: 



(79) 



In case of m\ 7^ 0: 



dY{ 



imY, 



1 vy> • 



(80) 
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sin(0i — 03) can be written in exponential form: 



sin 



"*03 _ e -i<t>l e i( t>3 



(81) 



Similarly one obtains: 



2l[ + 1 
2l[ + 3 



1/2 l[+l+h 



L 2 +l' 9 



x E ^ E ^ E (2) (2^i+l)(2L 2 +l) 1 / 2 

-^1 — |^1 + 1 — ^1 1 ^2 — 1^2~ ^| L'2 = \L2-\-l2 | 



x (24 + 1) 1/2 jc L2 (Z 2 , m' 2 ; / 2 , m 2 )C L2 (l, -1; Z 2 , m 2 )C^(L 2 , -m 2 - 1; l' 2 , m' 2 ) 



x 



m 



n2ir 
Y^-\e u ^Y^- 1 ^, 0i)e 2 ^ sinMflK 
2 



i + 2)(Zi - mi + l)] 1 / 2 C Ll (/'i + IX - 1; Zi,mO 



x 



'0 </0 

+ [(Z; + m' x + 1)(Z; + mi + 2)] 1 / 2 C Ll (Z' x + 1, mi + 1; Z 1; m x ) 

*»7T /*27T 




JO 



+ C i2 (l, 1; Z 2 , m 2 )C^ (La, -m 2 + 1; Z 2 , m' 2 ) 
x ([(Z; - mi + 2)(Z; - mi + 1)] 1/2 C Ll (Zi + 1, mi - 1; Zi, mi) 

Y^-\9i, 0i)r L f 2+1 (^i, 00 sin MM0i 

+ [(Zi + mi + l)(Zi + mi + 2)] 1 / 2 C Ll (Zi + 1, mi + 1; l u mi ) 

x /(JVi — 1, _/V" 2 + 1; wi, w 2 ; i/ + i/ - 2; L 2 )} (82) 

The expressions given here, although they look complicated, most of the terms vanish. They 
are very sensitive to errors. 



Especial case mi = 

The above expressions are general. For the case mi = 0, one may obtain simp 
sions using shorter recursion relations for the spherical harmonics Eq. (C.9) of 
can be easily obtained: 



er expres- 



29|. So it 
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iKEfiAV = (— l) M2 <5(Mi+M 2 , 0) 



(2Zi + l)(2Zi-l) 



{% + 1) 



-, 1/2 l[-l+h 

{2) (2L 1 +l) 1 / 2 C L ^l' 1 -l,m' 1 ;l 1 ,m 1 ) 

ii=|Zi-l-Zi| 



X 



(2l> 2 - 1) 



(Z 2 - m' 2 ){l> 2 + m' 2 )l 1/2 x ^ (2)5(Li) L , )(2L2 + 1)1/2 

L 2 =U'-l-i 2 | 



+ 



x C La (Z 2 - 1, m' 2 ; Z 2 , m 2 )I(N 1 -1,N 2 + 1; u u co 2 ; v + i/ - 2; L 2 ) 

1/2 ^+1+^2 

£ (2) (2^i + l) 1/2 (2L 2 + 1)^(^,4) 

L' 2 =|^+1_J 2 | 



(Z 2 -m 2 + l)(Z 2 + m 2 + l) 



(2l> 2 + 2) 

C L 2 (l 2 + 1, m 2 ; Z 2 , m 2 )I(N 1 -1,N 2 + 1; oo h oo 2 ; u + v' - 2; L' 2 ) } . 

m ' 1= 0, Zi = l,2 (83) 



and 



i' 9 +h 



I K E,e 1 ,2= 1 -(-l) Ml S(M 1 + M 2 ,0) {2) (2^2 + l) 1/2 ^ i2 (/ 2 ,m 2 ;Z 2 ,m 2 ) 

L2 = \l' 2 -h\ 
1/2 «i-2+«l 



X < (l[ - 1) 



(2/; + 1) 

(2Zi - 3) 



J2 (2) (2^i + l) 1//2 <K-£i, L 2 )C Ll (l[ - 2,m' 1 ;h,m 1 ) 

Li=|«i-2-h| 



+ /i E i2 \2L f 1 + l) 1 / 2 S(L' 1: L 2 )C L ^l' 1 ,m[-,h,m 1 ) 
ii=l'i-ii| 

x [/(iVx - 2, iV 2 ; Wl , w 2 ; i/ + i/; L 2 ) 3 ; 1, 1; L 2 , L 3 ) - /(A^, jV 2 ; Wl , w 2 ; i/ + i/ - 2; L 2 ) 
- /(TVi - 2, iV 2 + 2; w 1; w 2 ; i/ + v' - 2; L 2 )] . 

m[ = 0, Zi = l,2 (84) 

The derivative of the spherical harmonic with m' x = with respect to 4>i vanishes: 



Ike, foil) - 



(85) 



5.3 Construction of symmetry adapted configurations 



The configurations ss, pp, dd of S symmetry using complex Slater orbitals are constructed: 
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^ P = PP = Po(l)Po(2) - pi(l)p_i(2) - p_i(l)pi(2) = poPo - PiP-i - P-1P1 (86) 

the basis functions p\P-\ and p_iPi are energetically degenerated. A matrix element between 
two of these configurations has to be calculated as: 



VP 



H 



PP 



P0P0-P1P-1 -P-1P1 



H 



PoPo ~ PiP-i - P-1P1 



(87) 



The configuration dd of S symmetry is constructed: 



ipd — dd — dodo — dxd-i — d-\d\ + d 2 d_ 2 + d- 2 d 2 
the basis functions d\d_i and d~\d\ are degenerated, also d 2 d_ 2 and d_ 2 d 2 . 



(88) 



6. Calculations 



6.1 CI Calculations 

n 

In order to test the computer code the s,p, d- calculations of Table II of Weiss [5] were 
reproduced. The recalculated values are shown in Table 1. The CI codes using real and 
complex orbitals led to the same energy results. Both codes use quadruple precision, about 30 
decimal digits in our computer. The advantage of complex orbitals is that general algorithms 
can be developed and applied, as the angular integration, and general integral subroutines. 

We have tried first the full optimization of the CI wave function. Although the en- 
ergy results were very satisfactory, this technique had to be abbandoned because it become 
impracticable. 

In Table 2 we show a set of optimized exponents. They have been optimized in the 
field of about 10 configurations. Optimal exponents should be obtained by a generalized 
optimization procedure in the art made by Gaussian orbitals. We have used single and 
double, and half basis sets, these last were built up by multipliying the exponents of the 
single basis by 2 or dividing by 2. The double basis could be used until the exponent with 
value 20.0, further the quadruple precision was not enough. The half basis could be used 
only for some functions because linear dependence appeared (too similar exponents). 

The first 15 terms of the CI wave function in Table 3 are the same than the ones used 
by Weiss, they gave a good energy result. We also used the first 10 p-functions of Weiss. 
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Afterwards new terms are added systemstically. We show here the convergence of the 101- 
term CI wave function as an example. The convergence pattern of the CI wave function was 
very discouraging. In Table 4 we show our best results obtained using quadruple precision. 
They agree with the results of Sims and Hagstrom who used 124*real arithmetic and with 
the results of Bromley and Mitroy [34j using Laguerre functions. A 668 term CI calculation 
leads to —2.9027662664 a.u. which is about 1 milihartree accurate. The best calculation 
up to date and the largest expansion is the one of Sims and Hagstrom which is about 4 
microhartree above the exact value. 

Therefore we arrive to the same conclusion than Carrol, Silverstone and Metzger 32] 



and Kutzelnigg and Morgan [7J]: "For a two electron atom using conventional CI to achieve 
microhartree accuracy (they mentioned an accuracy of 20 microhartrees) requires enormous 
labor and sophisticacion, and indicates that further improvement in accuracy by brutal 
force CI would be extremely difficult, if not virtually impossible because of numerical linear 
dependence problems." 



6.2 Hy-CI Calculations 

We have taken the wave function of Weiss and have added the same functions containing 
ri2- The exponents were slighly optimized in a field of about 10 functions. The set of expo- 
nents of the ri2-configurations has been taken different than the purely CI configurations, 
see Table 1. They are in general slightly larger. In Table 5 a 62-term truncated Hy-CI wave 
function is shown. One can observe the faster convergence of the wave function leading to 
a result with 4 decimal digit accurate. Note that the computation time was 130 seconds, 
although the program code has not been refined to save time. 

Extensive calculations are shown in Table 6 using only s-, p- and ci-orbitals, i.e. L = 2, 
and quantum numbers N up to 16 allowed to obtain the energy —2.90372437683996324137 
a.u. with Virial factor 1.999999998787224, a result of nanohartree accuracy. 



Conclusions 



It has been shown that the CI wave function converges stremely slowly. The use of a 
set of optimized exponents improves the CI expansion. A shorter configuration expansion 
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is needed by the Hy-CI wave function to obtain a better result. Nanoharteee accuracy has 
been achieved using a set of s-, p- and rf-orbitals and r\2- This result can be improved by 
using a generalized Hy-CI basis set. 

Theoretically, the CI and Hy-CI calculations here presented agree with the ones of Sims 
and Hagstrom and this confirm the correctness of the Hamiltonian in Hylleraas coordinates 
and of the kinetic energy integral expressions. The results using the Hy-CI wave function are 
so promising, that it would be worth while the application of the method to larger systems. 
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Table 1: Repetition Weiss calculations [5 



Confs. E (a.u.) Virial 

Weiss 15 -2.8789 5525 8604 1.999 999 
Weiss 25 -2.9003 9002 0339 1.958 324 
Weiss 31 -2.9025 8300 1014 2.000 020 
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Table 2: Basis set used in the CI and in Hy-CI calculations. 



Orbital CI single CI double CI half Hy-CI single 



Is 


1.48 


2.96 0.74 


1.9 


2s 


1.48 


2.96 0.74 


1.9 


Is' 


3.7 




3.9 


2s' 


3.7 




3.4 


3s 


3.7 


7.4 1.85 


3.0 


4s 


4.7 


9.4 2.35 


5.4 


5s 


5.7 


11.4 2.85 


4.4 


6s 


6.7 


13.4 3.35 


5.5 


7s 


7.7 


15.4 3.85 


8.4 


8s 


8.7 


17.4 4.35 


8.8 


9s 


9.7 


19.4 4.85 


9.8 


10s 


10.7 


20.0 5.35 




lis 


11.7 


20.0 5.85 




12s 


12.7 


20.0 6.35 




13s 


13.7 


20.0 6.85 




14s 


14.7 


20.0 




15s 


15.7 


20.0 




2p 


2.7 


5.4 


3.0 


3p 


2.7 


5.4 


2.0 


2p' 


5.4 




5.0 


3p' 


5.4 




6.0 


4p 


3.3 


6.6 


3.8 


5p 


3.7 


7.4 


4.5 


6p 


4.3 


8.6 


5.2 


7p 


4.7 


9.4 


5.9 


8p 


5.3 


10.6 


6.5 


9p 


5.7 


11.4 


7.2 


Wp 


6.0 


12.0 


8.0 


lip 


6.3 


12.6 


8.7 


12p 


6.7 


14.0 


9.4 


13p 


7.0 


14.6 


10.2 


Up 


7.3 




11.0 


15p 


7.7 




11.6 


16p 


8.0 




12.4 


3d 


3.6 


7.2 


3.8 


id 


2.6 


5.2 


4.2 


5d 


4.7 


9.4 


4.6 


6d 


5.5 


11.0 


3.9 


Id 


6.3 


12.6 


4.1 


8d 


7.1 


14.2 


4.6 


9d 


7.8 


15.6 


5.7 


lOd 


8.6 


17.2 


6.4 


lid 


9.4 


18.8 


7.0 


12d 


10.0 


20.0 


7.6 


13d 


10.6 


21.2 


8.1 


14d 


11.4 




8.6 


15d 


12.2 




9.1 


16rf 


13.0 




9.6 
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Table 3: s-,p- and d- CI calculation on helium atom, 101 terms truncated wave function 



expansion. 



N Configuration 




Energy 




virial 




Difference 


1 


Isls 


-2 


8046 


0000 


0000 


2 


28041 






2 


ls2s 


-2 


8464 


5649 


2449 


2 


01265 


-0 


.041856492449 


3 


2s2s 


-2 


8528 


0522 


2663 


2 


03876 


-0 


.006348730214 


4 


Isls' 


-2 


.8749 


7028 


4519 


2 


00204 


-0 


.022165061856 


5 


2sls' 


-2 


.8753 


1452 


7730 


2 


00131 


-0 


.000344243211 


6 


WW 


-2 


.8762 


1971 


7811 


2 


00040 


-0 


.000905190081 


7 


ls2s' 


-2 


.8774 


6546 


4594 


2 


00105 


-0 


.001245746783 


8 


2s2s' 


-2 


8784 


4308 


9302 


2 


00006 


-0 


.000977624708 


9 


ls'2s' 


-2 


8784 


9065 


1589 


2 


00005 


-0 


.000047562287 


10 


2s'2s' 


-2 


8787 


4806 


7776 


1 


99986 


-0 


.000257416187 


11 


ls3s 


-2 


8787 


8660 


6389 


2 


00004 


-0 


.000038538613 


12 


2s3s 


-2 


8788 


2609 


4249 


2 


00009 


-0 


.000039487860 


13 


ls'3s 


-2 


8789 


0234 


0012 


2 


00003 


-0 


.000076245763 


14 


2s'3s 


-2 


.8789 


0276 


7328 


2 


00002 


-0 


.000000427316 


15 


3s3s 


-2 


.8789 


5525 


8604 


2 


00000 


-0 


.000052491276 


16 


ls4s 


-2 


.8789 


5642 


2107 


2 


00000 


-0 


.000001163503 


17 


2s4s 


-2 


.8789 


5663 


3396 


2 


00000 


-0 


.000000211289 


18 


3s4s 


-2 


.8789 


5751 


9018 


2 


00000 


-0 


.000000885622 


19 


4s4s 


-2 


8789 


7440 


1997 


1 


99999 


-0 


.000016882979 


20 


ls5s 


-2 


8789 


7575 


8637 


1 


99999 


-0 


.000001356640 


21 


2s5s 


-2 


8789 


7576 


2491 


1 


99999 


-0 


.000000003854 


22 


3s5s 


-2 


8789 


7720 


3435 


1 


99999 


-0 


.000001440944 


23 


4s5s 


-2 


8789 


8084 


9335 


1 


99999 


-0 


.000003645900 


24 


5s5s 


-2 


8789 


8702 


4993 


1 


99998 


-0 


.000006175658 


25 


ls6s 


-2 


8789 


8769 


6621 


1 


99998 


-0 


.000000671628 


26 


2s6s 


-2 


8789 


8772 


2861 


1 


99998 


-0 


.000000026240 


27 


3s6s 


-2 


.8789 


8826 


1084 


1 


99998 


-0 


.000000538223 


28 


4s6s 


-2 


.8789 


8868 


8728 


1 


99998 


-0 


.000000427644 


29 


5s6s 


-2 


8789 


8975 


9206 


1 


99998 


-0 


.000001070478 


30 


6s6s 


-2 


8789 


9115 


9605 


1 


99998 


-0 


.000001400399 


31 


ls7s 


-2 


8789 


9262 


2222 


1 


99999 


-0 


.000001462617 


32 


2s7s 


-2 


8789 


9262 


9635 


1 


99999 


-0 


.000000007413 


33 


3s7s 


-2 


8789 


9327 


4255 


1 


99998 


-0 


.000000644620 


34 


4s7s 


-2 


8789 


9335 


5895 


1 


99999 


-0 


.000000081640 


35 


5s7s 


-2 


8789 


9346 


3799 


1 


99998 


-0 


.000000107904 


36 


6s7s 


-2 


.8789 


9377 


7455 


1 


99998 


-0 


.000000313656 


37 


7s7s 


-2 


.8789 


9419 


5042 


1 


99998 


-0 


.000000417587 


38 


ls8s 


-2 


.8789 


9565 


2668 


1 


99999 


-0 


.000001457626 


39 


2s8s 


-2 


.8789 


9565 


4581 


1 


99999 


-0 


.000000001913 


40 


3s8s 


_2 


8789 


9600 


1814 


1 


99999 


_Q 


UVI'JV'UVI.).|Ui.).) 


41 


4s8s 


-2 


.8789 


9606 


3092 


1 


99999 


-0 


.000000058278 


42 


5s8s 


-2 


.8789 


9608 


4638 


1 


99999 


-0 


.000000021546 


43 


6s8s 


-2 


8789 


9611 


7732 


1 


99999 


-0 


.000000033094 


44 


7s8s 


-2 


8789 


9624 


1954 


1 


99999 


-0 


.000000124222 


45 


8s8s 


-2 


8789 


9639 


3931 


1 


99999 


-0 


.000000151977 


46 


2p2p 


-2 


8980 


1066 


3847 


1 


.99720 


-0 


.019014269916 


47 


2pSp 


-2 


8983 


0353 


5068 


1 


99927 


-0 


.000292871213 


48 


3p3p 


-2 


8999 


8445 


0188 


2 


00053 


-0 


.001680915128 



27 



Continuation Table 3. 



N Configuration 




Energy 




virial 


Difference 


49 


2p2p' 


-2 


.9000 


1939 


5255 


2.00040 


-0.000034945067 


50 


3p2p' 


-2 


.9000 


2572 


8711 


2.00030 


-0.000006333456 


51 


2p'2p' 


-2 


.9003 


0029 


7330 


2.00012 


-0.000274568619 


52 


2p3p' 


-2 


.9003 


0993 


6342 


2.00014 


-0.000009639012 


53 


3p3p' 


-2 


.9003 


4614 


9487 


2.00001 


-0.000036213145 


54 


2p'3p' 


-2 


.9003 


4931 


2787 


2.00002 


-0.000003163300 


55 


3p'3p' 


-2 


.9004 


2459 


8409 


2.00001 


-0.000075285622 


56 


2p4p 


-2 


.9004 


3166 


9931 


2.00003 


-0.000007071522 


57 


3p4p 


-2 


.9004 


3221 


0413 


2.00004 


-0.000000540482 


58 


4p4p 


-2 


.9004 


3399 


3582 


2.00004 


-0.000001783169 


59 


2p5p 


-2 


.9004 


3400 


6808 


2.00004 


-0.000000013226 


60 


3p5p 


-2 


.9004 


3508 


9064 


2.00005 


-0.000001082256 


61 


4p5p 


-2 


.9004 


4229 


5990 


2.00002 


-0.000007206926 


62 


5p5p 


-2 


.9004 


4249 


5331 


2.00002 


-0.000000199341 


63 


2p6p 


-2 


.9004 


4266 


5968 


2.00002 


-0.000000170637 


64 


3p6p 


-2 


.9004 


4268 


0235 


2.00002 


-0.000000014267 


65 


4p6p 


-2 


.9004 


4712 


1025 


2.00002 


-0.000004440790 


66 


5p6p 


-2 


.9004 


5007 


1856 


2.00003 


-0.000002950831 


67 


6p6p 


-2 


.9004 


5688 


7282 


2.00002 


-0.000006815426 


68 


2p7p 


-2 


.9004 


5688 


7957 


2.00002 


-0.000000000675 


69 


3p7p 


-2 


.9004 


5691 


7483 


2.00002 


-0.000000029526 


70 


4p7p 


-2 


.9004 


6603 


9818 


2.00001 


-0.000009122335 


71 


5p7p 


-2 


.9004 


6617 


5403 


2.00001 


-0.000000135585 


72 


6p7p 


-2 


.9004 


6618 


6237 


2.00001 


-0.000000010834 


73 


7p7p 


-2 


.9004 


7220 


2637 


2.00001 


-0.000006016400 


74 


2p8p 


-2 


.9004 


7221 


0819 


2.00001 


-0.000000008182 


75 


3p8p 


-2 


.9004 


7221 


5238 


2.00001 


-0.000000004419 


76 


4p8p 


-2 


.9004 


7250 


5447 


2.00001 


-0.000000290209 


77 


5p8p 


-2 


.9004 


7255 


4026 


2.00001 


-0.000000048579 


78 


6p8p 


-2 


.9004 


7428 


5314 


2.00001 


-0.000001731288 


79 


7p8p 


-2 


.9004 


7683 


0335 


2.00001 


-0.000002545021 


80 


8p8p 


-2 


.9004 


7916 


7253 


2.00001 


-0.000002336918 


81 


3d3d 


-2 


.9020 


6435 


2165 


1.99950 


-0.001585184912 


82 


3d4d 


-2 


.9021 


5865 


5625 


1.99974 


-0.000094303460 


83 


4d4d 


-2 


.9025 


6734 


4690 


1.99988 


-0.000408689065 


84 


3d5d 


-2 


.9025 


7143 


9942 


1.99990 


-0.000004095252 


85 


4d5d 


-2 


.9025 


7753 


1336 


1.99995 


-0.000006091394 


86 


5d5d 


-2 


.9026 


6778 


7092 


2.00001 


-0.000090255756 


87 


3d6d 


-2 


.9026 


7329 


8893 


2.00002 


-0.000005511801 


88 


4d6d 


-2 


.9026 


7446 


9528 


2.00002 


-0.000001170635 


89 


5d6d 


-2 


.9026 


7512 


9159 


2.00001 


-0.000000659631 


90 


6d6d 


-2 


.9026 


9774 


0952 


2.00002 


-0.000022611793 


91 


3d7d 


-2 


.9026 


9775 


3228 


2.00002 


-0.000000012276 


92 


idld 


-2 


.9027 


0076 


2033 


2.00003 


-0.000003008805 


93 


5d7d 


-2 


.9027 


0111 


8387 


2.00003 


-0.000000356354 


94 


6dJd 


-2 


.9027 


0375 


0417 


2.00002 


-0.000002632030 


95 


Idld 




.9027 


1095 


2280 


2.00002 


-0.000007201863 
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Continuation Table 3, 



N Configuration 


Energy virial 


Difference 


96 


3dSd 


9 QD97 1 1 nfi 4338 9 nnnn9 




97 


4d&d 


-2.9027 1112 0091 2.00002 


-0.000000055753 


98 


5dSd 


-2.9027 1223 1094 2.00002 


-0.000001111003 


99 


6d8d 


-2.9027 1244 9861 2.00002 


-0.000000218767 


100 


TdU 


-2.9027 1499 4392 2.00002 


-0.000002544531 


101 


Sd&d 


-2.9027 1789 6028 2.00002 


-0.000002901636 


Exact 




-2.9037 2443 1992 


-0.001006535964 



Calculation time 537 seconds. 



Table 4: Comparison of extensive spd-CI calculations on helium atom. 



Best CI 



'max 


max 


Confs. 


This work 


■^Y max 


Confs. Sims and Hagstrom f3] iV or b 


Bromley and Mitroy 





i.". 


296 


-2.879 028 6809 


21 


470 


-2.879 028 7565 


44 


-2.879 028 760 


l 


10 


497 


-2.900 515 8544 


21 


854 


-2.900 516 2199 


80 


-2.900 516 228 


2 


10 


668 


-2.902 766 2664 


21 


1221 


-2.902 766 8053 


115 


-2.902 766 823 


18 


i.". 








4699 


-2.903 720 091 9 


465 


-2.903 712 786 



[34] 



In this work were used the exponents of Table 1. Sims and Hagstrom used 4.10 and 25.0 for 

s-orbitals, 

3.05 and 40.5 for p-orbitals and 3.50 and 40.5 for ci-orbitals. 
Bromley and Mitroy used Laguerre functions with the exponent 8.6 for s-orbitals, 11.6 for 
p-orbitals and 14.4 for ci-orbitals. -/V orD is the number of orbitals. 
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Table 5: s-,p- and d- Hy-CI calculation on helium atom with optimization of the orbital 

exponents from Table 1. 



N Configuration 






Energy 








Difference 


1 


lsls 


_2 


8040 


ooi id 


001 1() 


( )( it i( ) 


( )( )()( ) 


9 980405 




2 


ls2s 


_2 


s \{\ ; 


5049 


2448 


7880 


8848 


9 01 9655 


fl 041 8564Q944Q 
u . ut: x oou^ty ^TiTiy 


3 


2s 2s 


_2 


8528 


0o22 


2002 


7000 


0893 


9 0^8758 


00ft'}487 i; l0914 
-U.UUOO^tO / OUZ 14 


4 


lsls' 


_2 


8749 


7028 


4ol8 


0921 


8031 


9 00904'} 

Z.UUZU4:0 


0991 fi^Ofil 85ft 

-U.UZZ1UOUOXOOO 


c 



ZS1S 


_2 


.8753 


1 A 

140Z 


7740 


4383 


9933 


9 001 30fi 

ii.UUXU uu 


fl 00034494391 9 

U . UUUOt:t:Z-XOZ X Z 


6 


is is 


_2 


.8762 


1971 


7810 


5121 


5735 


9 000400 

Zi.UU U1UU 


fl 0009051 Q0080 


7 


Is 2s' 


_2 


8774 


0o40 


4o94 


0870 


3077 


9 001 049 
z.uuiU4y 


fl 001 94^741^784 
-U.UU1Z410 / 40 / 04 


8 


2s2s' 


_2 


8784 


4308 


9302 


1401 


8945 


9 OOOOft^i 


000077ft94708 
-u.uuuy 1 f uz4 / uo 


9 


ls^s' 


_2 


.8784 


9005 


L>88 


5673 


0022 


9 000045 
Zi.UU uuu 


000047569986 


10 


2s'2s / 


_2 


8787 


4800 


7775 


0347 


■>228 


1 QQOS59 


fl 00095741 ftl 87 
-\.).UVJy.)Z,0 1 41010 / 


i i 


ls3s 


_2 


.8787 






3280 


6871 


9 nonn^fi 

a.UUUUuU 


00003853861 4 


12 


2s3s 


_2 


8788 


2009 


4249 


0916 


2944 


9 000087 


fl 000030487860 
u.uuuuoyto 1 ouu 


13 


ls'3s 


_2 


8789 


0234 


0012 


1737 


8000 


9 000095 
z .uuuu zo 


fl 000076945763 
-U.UUUU 1 OZ40 1 OO 


14 


2s'3s 


_2 


8789 


0270 


7327 


8845 


3030 


9 000099 
z.uu uuzz 


00000049731 ft 
-U.UUUUUU4Z ( OlO 


10 


OSoS 


_2 


.878!) 


OOZO 


OOUo 


6072 


3426 


1 qqqqqq 


fl 0000594Q1 97ft 
u .uuuuu Zi^ty xz 1 u 


l(i 


lslsri2 


_2 


.<)030 


3855 


5047 


6765 


9888 


9 000060 


0940839QR444 


17 


1s2sT"12 


_2 


9030 


3903 


0490 


7037 


1147 


1 9Q9Q05 
x.yyyyuo 


000001 07544 3 
-U.UUUUU1U 1 0440 


1 Q 
10 


2s2sT"i2 


_2 


.9032 




ZOOO 


4800 


4038 


9 0000^0 


00091 898906ft 
~y.>.\jyjy.iz, xoyoz<uuu 


1!) 


lslsVi2 


_2 


.9033 


8299 


0275 


5132 


6846 


9 000091 
a.uu uuy x 


0001 9438371 9 


20 


As Is ri2 


_2 


9034 


1178 


0230 


9870 


1758 


9 00001 9 

Z.UU UUXZ 


00009R7R3Q61 
-u.uuuuzo / ooyox 


21 


ls'lsVi2 


_2 


.9034 


3722 


4910 


8831 


4415 


1 9Q9Q89 
x .yyyyoy 


000095444704 

U . UUUUZOtitit: 


22 


Is2s'ri2 


_2 


9034 


3950 


8142 


1374 


1002 


x .yyyy 1 z 


fl 000009983901 
-u.uuuuuzzoozux 


23 


2s2sVi2 


-2 


9034 


0815 


( i 198 


5144 


023'-> 


x .yyy you 


000098ft4935ft 


24 


ls'2sVl2 


_2 


9034 


0820 


5099 


9014 


1845 


1 QQQG^ft 
x .yyyyou 


n nonoono r i4ftni 

U . UUUUUUU04UU X 


OK 
ZO 


ZS Z,S Tj2 


_2 


.9034 


7Q1 1 
( oil 


A KM 
40 01 


4096 


4381 


1 9Q9Q59 
x .yyyyou 


000004909451 


20 


Is3sri2 


-2 


.9034 


7048 


8273 


0412 


1726 


1 QQQQQ9 
x .yyyyyz 


00000T}7'}799 
u .yyjyjyjyjoo 1 1 zz 


27 


2s3sri2 


_2 


9034 


7957 


2554 


5670 


0172 


9 000009 
z.uuuuuz 


00000^084989 
-U.UUUUUOU04ZOZ 


ZO 


ls'3sri2 


_2 


.9034 


oloo 


O ( Zo 


1881 


1971 


9 00000'} 

a.UU uuuo 


000001 7841 69 

U.UUUUUl 1 041UJ 


29 


2s'3.sri2 


-2 


.9034 


8172 


3708 


1817 


9421 


9 000009 

i.UU UUUi 


000000 '}6 704 5 

U.UUUUUUOU 1 U4:0 


Qfl 
OU 


3s3sr^2 


-2 


.9034 


oiyi 


UZZo 


3344 


9701 


9 00000'^ 
Zi.UU uuuo 


fl 0000001 86460 

U.UUUUUUXOU4:UU 


31 


Ap2p 


_2 


.9030 


9490 


9181 


3235 


6472 


1 9Q9Q40 


00091 9QQ8Q53 
-y.>.ijyjy.>^ xzyyooot> 


32 


2p3p 


-2 


9037 


0333 


9014 


5345 


8498 


1 0Q0Q7^t 

1 .yyyy 1 


0000084998^^ 
-u.uuuuuo4zyooo 


33 


3p3p 


-2 


.9037 


1382 


9571 


8228 


7221 


1.999986 


-0.000010490557 


34 


2p2p' 


-2 


.9037 


1407 


2375 


9450 


1950 


1.999985 


-0.000000242804 


35 


3p2p' 


-2 


.9037 


1429 


5379 


2885 


5478 


1.999989 


-0.000000223003 


36 


2p'2p' 


-2 


.9037 


1572 


6120 


9390 


5519 


1.999990 


-0.000001430742 


37 


2p3p' 


-2 


.9037 


1588 


5177 


7166 


9859 


1.999993 


-0.000000159057 


38 


3p3p' 


-2 


.9037 


1634 


2690 


0297 


1140 


1.999994 


-0.000000457512 


39 


2p'3p' 


-2 


.9037 


1644 


9779 


9809 


6406 


1.999993 


-0.000000107090 


40 


3p'3p' 




.9037 


1644 


9785 


9439 




1.999993 


-0.000000000006 
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Continuation Tabic 5. 



N Configuration 



Energy 



virial 



Difference 



41 2p2pr 12 

42 2p3pr 12 

43 3p3pri2 

44 2p2pVi2 

45 3p2p'n 2 

46 2p'2p'r 12 

47 2p3pVi2 

48 3p3p'ri2 

49 2p'3pVi 2 

50 3p'3pVi2 

51 3d3d 

52 MAd 

53 4d4d 

54 Mbd 

55 4d5d 

56 5d5d 

57 3d3dri2 

58 3d4dri2 

59 AdAdr 12 

60 3d5dri2 

61 4d5dri2 

62 5d5dri2 
Exact 



-2.9037 2173 2417 9137 9221 1.999989 
-2.9037 2182 0699 2263 4691 1.999990 
-2.9037 2214 9854 2463 6331 1.999993 
-2.9037 2227 8371 0330 7986 1.999994 
-2.9037 2237 5963 4585 9805 1.999994 
-2.9037 2265 5745 0023 6384 1.999994 
-2.9037 2272 5043 9425 8433 1.999994 
-2.9037 2278 0206 3073 0904 1.999995 
-2.9037 2287 1724 9868 5604 1.999996 
-2.9037 2311 0207 5489 3605 1.999997 
-2.9037 2311 5896 1644 4399 1.999997 
-2.9037 2325 3694 3809 8700 1.999996 
-2.9037 2341 0977 7800 8564 1.999997 
-2.9037 2363 2223 7418 5292 1.999998 
-2.9037 2368 3956 6486 0703 1.999998 
-2.9037 2371 6568 3659 1072 1.999999 
-2.9037 2375 1435 0361 1170 1.999999 
-2.9037 2375 1898 7744 7923 1.999999 
-2.9037 2376 3066 1777 0715 1.999999 
-2.9037 2377 6768 3902 6301 1.999999 
-2.9037 2377 8464 9764 7830 1.999999 
-2.9037 2378 4817 0105 8051 1.999999 
-2.9037 2437 7034 1195 9831 



-0.000005282632 
-0.000000088281 
-0.000000329155 
-0.000000128517 
-0.000000097592 
-0.000000279782 
-0.000000069299 
-0.000000055162 
-0.000000091519 
-0.000000238483 
-0.000000005689 
-0.000000137798 
-0.000000157283 
-0.000000221246 
-0.000000051733 
-0.000000032612 
-0.000000034867 
-0.000000000464 
-0.000000011167 
-0.000000013702 
-0.000000001697 
-0.000000006352 
-0.000000592217 



Calculation time: 130 seconds. 0.5 microhartree accuracy. 
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Table 6: Comparison of extensive Hy-CI calculations on helium atom using single basis set. 



-^max 


N 

1 v max 


Confs. This work 


Accuracy 





15 


378 


-2.9034 9813 4711 9664 5493 


milihartree 


1 


16 


482 


-2.9037 2426 1717 2135 6263 


microhartree 


2 


16 


692 


-2.9037 2437 6839 9632 4137 


nanohartree 


Exact 






-2.9037 2437 7034 1195 9831 12 





These values are obtained in QP. In this work the exponents of Table 1 were used. 
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